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Regular factorial designs with randomization restrictions are widely 
used in practice. This paper provides a unified approach to the con- 
struction of such designs using randomization defining contrast sub- 
spaces for the representation of randomization restrictions. We use 
finite projective geometry to determine the existence of designs with 
the required structure and develop a systematic approach for their 
construction. An attractive feature is that commonly used factorial 
designs with randomization restrictions are special cases of this gen- 
eral representation. Issues related to the use of these designs for par- 
ticular factorial experiments are also addressed. 

1. Introduction. The use of factorial experiments in situations which re- 
quire randomization restrictions (e.g., block designs or split-plot designs) 
has been illustrated in the literature for many years [see, e.g., Cochran and 
Cox (1957), Chapters 6-8 and Addelman (1964)]. More recently, the con- 
struction of designs under various optimality and efficiency criteria has been 
discussed for blocked factorial and fractional factorial designs by Bisgaard 
(1994), Sitter, Chen and Feder (1997), Sun, Wu and Chen (1997), Chen and 
Cheng (1999) and Cheng, Li and Ye (2004); for factorial and fractional fac- 
torial split-plot designs by Huang, Chen and Voelkel (1998), Bingham and 
Sitter (1999) and Stapleton, Lewis and Dean (2009); for strip-plot designs 
by Miller (1997); and for split-lot designs by Mee and Bates (1998) and 
Butler (2004). Although these various designs maintain the same factorial 
treatment structure, their randomization structures are different. 

Bingham et al. (2008) proposed the use of randomization defining contrast 
subgroups (RDCSGs) to describe the randomization structure of factorial 
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designs. These can be viewed as generalizations of block defining contrast 
subgroups [see, e.g., Sun, Wu and Chen (1997)]. In this article, we propose 
a projective geometric generalization of the RDCSGs, which we refer to as 
randomization defining contrast subspaces (RDCSSs). 

For complicated randomization structures, the existence of desirable de- 
signs can be difficult to determine; for example, Bingham et al. (2008) were 
forced to search over the set of all possible designs to find those satisfying the 
required restrictions. Here, we develop theoretical results for the existence 
of factorial designs with given randomization restrictions within a unified 
framework, and we provide a direct method of construction of such designs. 
The theory for the existence and construction of such designs is given in 
terms of full factorial designs, but we show in Section 6 that the results 
apply equally well to fractional factorial designs. 

In the next section, a brief description of RDCSSs is presented. Section 
3 considers the impact of the randomization restrictions on data analysis 
and motivates the preference for nonoverlapping RDCSSs over the various 
stages of randomization. In Section 4, theoretical results for the existence 
of designs with randomization restrictions are developed. The theoretical 
framework is constructive, thereby allowing for the identification of designs 
in practical settings. Methods of construction of such designs for specific 
factorial experiments are developed in Section 5. The article finishes with a 
discussion of fractional factorial designs in Section 6 and concluding remarks 
in Section 7. 

2. Randomization-defining contrast subspaces. A natural design choice 
for a p-factor experiment is a completely randomized design with trials per- 
formed in a random order. However, this is not always feasible as, for exam- 
ple, when the experimental units are not homogeneous (as in block designs), 
the levels of some factors cannot be changed as rapidly as that of others 
(as in split-plot designs) or when subsets of factors must be held fixed at 
different stages of experimentation (as in split-lot designs). 

Several different approaches have been taken in the construction of fac- 
torial designs with randomization restrictions. For example, blocking in fac- 
torial and fractional factorial experiments have been studied using Abelian 
groups and vector spaces [see, e.g., Dean and John (1975), Bailey (1977) and 
Voss and Dean (1987)], and using finite geometries [see, e.g., Bose (1947), Sri- 
vastava (1987) and Mukerjee and Wu (1999)]; the construction of strip-plot 
designs via latin square fractions has been discussed by Miller (1997); split- 
plot designs for factorial treatments have been constructed using computer 
search and group theory by Addelman (1964), Bingham and Sitter (1999), 
Bisgaard (2000), Huang, Chen and Voelkel (1998) and Stapleton, Lewis and 
Dean (2009). For the construction of split-lot designs, linear graphs were 
used by Taguchi (1987), combinations of linear graphs and cyclic groups by 
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Mee and Bates (1998) and a grid-representation technique by Butler (2004). 
In this paper, we use finite projective geometry to study q p factorial designs 
with randomization restrictions for prime or prime power q. 

As in Dey and Mukerjee (1999), Chapter 8, let b be a p-dimensional 
pencil over the Galois field GF(q). For a(^ 0) G GF(q), b and ab represent 
the same pencil carrying q — 1 degrees of freedom. A pencil b represents an 
r-factor interaction if b has exactly r nonzero elements. The set of all p- 
dimensional pencils over GF(q) forms a (p— 1) -dimensional finite projective 
geometry and is denoted by PG{p — l,q). There are {q p — l)/{q — 1) points 
in PG(p — l,q), where points correspond to pencils. 

Since two-level factorial designs are the most common in practice, we 
focus on q = 2, although most of the results in this article hold for a prime 
or prime power q. For q = 2, a pencil 6 with r nonzero elements corresponds 
to a unique r-factor interaction in a 2 P factorial design with a single degree 
of freedom. Thus, the set of all effects (excluding the grand mean) of a 2 P 
factorial design is equivalent to the set of points in PG(p — 1,2) and will be 
referred to as the effect space V. 

The restrictions on the randomization of experimental runs is equivalent 
to the grouping of experimental units into sets of trials. We consider the 
usual approach of using independent effects from V to define the groupings. 
Blocked factorial designs, for example, use 2* (t < p) combinations of t inde- 
pendent effects from V to divide 2 P treatment combinations into 2* blocks. 
These factorial effects are then completely confounded with block effects 
and represent t randomization restriction "factors." The set S of all nonnull 
linear combinations of these t randomization restriction factors in V, over 
GF(2), forms a (t — 1) -dimensional projective subspace of V = PG(p — 1,2), 
which we call a randomization defining contrast subspace (RDCSS). 

Example 1. Consider an experiment arranged as a 2 5 factorial split-plot 
design, where A, B are whole-plot factors and C, D, E are sub-plot factors. 
The effect space is V = (A,B,C,D,E) and the RDCSS that imposes the 
randomization restrictions is S = {A,B) = {A,B,AB}, where the notation 
(ai, . . . , ak) denotes the projective space spanned by cti, . . . , a& (or, for two- 
level factorial designs, the set of all interactions of a±, . . . ,%). The notation 
A, AB, and so on, represent the main effect of A, the interaction of A 
and B, and so on. Since there are t = 2 whole-plot factors, the set of all 
experimental units is partitioned into 2* = 4 subsets (batches, whole-plots, 
etc.), and each subset consists of 2 5 /2 2 = 8 experimental units. The four 
subsets, say B\, £>2, B3 and #4, consist of experimental units corresponding 
to (9A(i),0B(i)) = (0, 0), (0, 1), (1,0) and (1,1), respectively, where, 9$(i) is 
the ith row of the column corresponding to the factorial effect 5 in the model 
matrix X [see, e.g., Bingham et al. (2008) and page 26 in Ranjan (2007)]. 
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There may be more than one stage of randomization restriction in a fac- 
torial experiment where the randomization structure can be characterized 
by its RDCSSs. For a 2 P factorial design with m stages of randomization, 
the m RDCSSs can be denoted by the projective subspaces Si, ... , S m con- 
tained in the effect space V. For each i = 1, . . . ,m, the size of Si is 2** — 1 
with < ti < p. Then, at stage i, the experimental units are partitioned into 
| Si | + 1 sets (e.g., batches, whole-plots and blocks) due to Si, where the size 
of each set is (\V\ + l)/(|Si| + 1) = 2 p ~ ti . 

Example 2. Consider a 2 5 factorial experiment with randomization 
structure defined by a strip-plot design [Miller (1997)], where the row config- 
urations are defined by a 2 2 design in factors A, B and the column configura- 
tions by a 2 3 design in factors C, D, E. In this setting, the two RDCSSs are 
Si = {A, B) and S 2 = {C,D,E), and the effect space is V = {A, B, C, D, E) . 

Although the treatment structure for Examples 1 and 2 are the same, the 
randomization restrictions induce different error structures, thereby having 
an impact on the analysis [see, e.g., Milliken and Johnson (1984), Chapter 

4]- 

3. Modeling issues. In this section, we show how the distribution of the 
least squares estimator of a factorial effect is related to the RDCSSs. This 
motivates a design strategy and, in particular, suggests a preference for 
nonoverlapping RDCSSs. 

3.1. Model. Consider a single-replicate 2 P factorial design with the linear 
regression model as the response model of interest; that is, 

(1) Y = X(3 + e, 

where X denotes the n x 2 P model matrix and (3 = (/3q, 0\, . . . , fop-i)' is the 
2 P x 1 vector of parameters corresponding to the grand mean, the factorial 
main effects and interactions. In this section, the factor levels and 1 are 
recoded as +1 and —1, respectively. For a single-replicate 2 V factorial experi- 
ment, n = 2 P , and the model matrix, X , is a Hadamard matrix which satisfies 
X'X = nl n , where I n is an n x n identity matrix. Without loss of generality, 
any p independent columns of X can be selected to represent the main effect 
contrasts of the p factors [see, e.g., Dean and Voss (1999), Section 15.6]. On 
rearranging the columns of X, let X be {co, c\, . . . , Cp, c p +i, . . . , c n _i}, where 
the column vector, cq, consists of all l's corresponding to the grand mean, 
columns c\ , . . . , c p refer to the independent main effect contrasts of the p 
factors and the remaining columns, c p +i, . . . ,c n _i, represent the interaction 
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contrasts obtained as element- wise products of subsets of c±, . . . ,c p . The vec- 
tor Y denotes the vector of response variables and e the error vector, whose 
distribution is discussed below. 

For a factorial design with m stages of randomization and RDCSSs de- 
noted by Si, i = 1, . . . , m, the error vector e can be written as a sum of m + 1 
independent error vectors, e = e$ + e\ + • • • + e m . The n X 1 vector eo denotes 
the replication error vector, and £j (1 < i < m) is the error vector associated 
with the randomization restriction characterized by Si. Since Si creates a 
partition of the n experimental units into 2 li batches (blocks, whole-plots, 
etc.), the error vector £j can be written as iVjej, where ej is a 2 f ' x 1 vector 
of errors associated with each of the 2** batches; the n x 2** matrix Ni is 
called the ith incidence matrix and has elements defined as follows: 

(Ni)rl = 1) if the rth experimental unit belongs to the /th batch 

(2) at the ith stage of randomization, 

= 0, otherwise, 
for r = 1, . . . , n and I = 1, . . . , 2 * . Thus, e can be written as 

(3) e = fr + £iH h£ m = eo + iVieiH h N m e m , 

where eo, ex, . . . , e m are independent, and eo = £q ~ N(0 n , I n a 2 ), where n is 
annxl vector of zeros and I n is an n x n identity matrix. For i = 1, . . . ,m, 
we assume that ej ~ iV(0 2 i, , I 2 tiCrf), where 2 ti is a vector of 2* 1 zero elements. 
It follows that e ~ N(0 n , T, y ), where 

m 

(4) E 1/ = < 7 2 7 n + 5^<7 i 2 ^. 

i=l 

Lemma 1 gives some properties of the incidence matrices which are needed 
for subsequent results. The proof of Lemma 1 is straightforward and omitted. 

Lemma 1. Consider a 2 P factorial design with n = 2 p runs and m levels 
of randomization restrictions defined by RDCSSs Si, ... , S m . Let N\, . . . , N m 
be the incidence matrices corresponding to Si's as defined in (2). Then, 
N-Ni = nil 2 ti , where m = 2 p ~ li is the number of 1 's in each column of Ni, 
and I 2 t l is a 2 li x 2 li identity matrix. 

3.2. Distribution of effect estimators. The most natural way to estimate 
the factorial effect parameters is to use the generalized least squares (GLS) 
estimator (3 = (X'Y,~ 1 X)~ 1 X'Y,~ l Y . Due to the assumptions on the error 
vectors, (3 ~ N(f3, Var(/3)). For a single-replicate 2 P factorial design, the GLS 
estimator coincides with the ordinary least squares (OLS) estimator (3 = 
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(X'X) 1 X'Y. This follows since X is an n x n matrix and, consequently, 
X'X = nl implies that X' = nX~ l . Then, (3 = n~ l X'Y, and 

$ = (X'E^X^X'E^Y = X'^yiX'y^X'^Y = n~ l X'Y = (j. 

Thus, the variance-covariance matrix of the factorial effect estimator (3 is 
Var(/3) = Var(/3) = {X'^ y X)/n 2 , and 

(5) P~N(/3,X%X/n 2 ), 

where T, y is defined in (4). 

If SVs are projective subspaces in V , it is sometimes possible to select the 
RDCSSs Si, ... , S m so that Sij = Si n Sj = 4> for all i ^ j. It will become 
clear later that these cases are of specific interest to practitioners since, in 
this case, all the effects in Si have variances that are functions of a 2 only. 
However, when this condition does not hold (see Theorem 2, below), the 
variance of the effects in Sij are impacted by both of and a 2 . 

Theorem 1. In a single-replicate 2 P factorial design, let the random- 
ization restrictions be defined by Si,...,S m contained in V. Then, for any 
two effects Ej and E k in the effect space V , the corresponding parameter 
estimators Pej and (3E k have independent normal distributions. 

Proof. Since (3 has a multivariate normal distribution, it is enough to 
show that cov(Pej , PE k ) = 0. Let Ej denote the factorial effect corresponding 
to jth column Cj (j > 1) of X, and let Si , . . . , S m be subspaces in V. If column 
Cj of X is orthogonal to the randomization restrictions in Si, then c'jNi = 0. 
Whereas, for each i = 1, . . . ,m, if Ej £ Si, then, as in Dean (1978), Cj may 

be expressed as Cj = Xia^ , where afj is a 2*' x 1 vector, and aj' 1 = with 
1 a vector of l's. Also, since c'jCk = for all j ^ k, it follows from Lemma 1 

that afp ap = for all j ^ k and i = 1, . . . , m. 

Define the index set Tj^ = {i : 1 < i < m, Ej S Si and E^ G Si}. Then, us- 
ing the equivalence of the GLS and OLS estimators for (3 together with (4), 
(5) and Lemma 1, we have 



in 



= £ <%af 'NlNiNlNia.% 



E2 2 («)' (i) n 



and the proof follows. □ 



Next, the main result of this section establishes the link between the vari- 
ances of the effect estimators for a 2 P full factorial design and the RDCSSs. 
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Theorem 2. In a single-replicate 2 P factorial design, let the randomiza- 
tion restrictions be defined by Si,..., S m in V . For an effect E E V , define 
an index set {Te, E E V} such that Te = {i : 1 < i < m, E G Si} . Then, 

f f! + £ tfE e{S 1 u---uS m }, 
Var(/3 E ) = I U 2 ieT E n 

— , pe?\{S 1 u-uS ra } ) 

v n 

where a 2 and of are the variances corresponding to 6q and Ei, respectively. 

Proof. Let E be the effect corresponding to column c in X. Following 
arguments similar to those of the proof of Theorem 1, if E E Si, we have 

(6) c'c = a^'NlNia^ = n;a (i) 'a« = n 

and, if E Si, then c'Ni = 0. Thus, 

n 2 Var E ) =ncr 2 +^ oja^' ' NtNiNtNiO® 

= no 2 + a\n ■ rii [using Lemma 1 and (6)] 

as required. On the other hand, if E G V \ {Si U • • • U S m }, Te is the empty 
set and Var {(3e) = a 2 /n. □ 

We note that the results in Theorems 1 and 2 can be extended to regular 
designs with factors having q levels for prime or prime power q. This is done 
by including in the model matrix a set of orthogonal contrasts for each effect, 
with each contrast scaled so that c'c = n. 

A common strategy for the analysis of unreplicated factorial designs is the 
use of half-normal plots [Daniel (1959)]. These require that effects whose es- 
timates appear on the same plot must be independent and have the same 
variance. By Theorems 1 and 2, when Sij = Si n Sj = <p,\/i ^ j, m sepa- 
rate half-normal plots (plus an additional plot if V \ {UilLi <Si} 7^ <f>) can De 
constructed to assess the significance of the effects. On the other hand, if 
Sij ^ (f> for some then the effects in Sij will have variances that are 
linear combinations of o~f and a 2 , creating a larger number of smaller sets 
of effects having the same variance. In this case, assessment of the signif- 
icance of effects in Sij may have to be sacrificed due to lack of degrees of 
freedom [Schoen (1999)]. Thus, designs with disjoint RDCSSs (Sij = 4> for 
all i ^ j) are preferred. From a practical standpoint, determining whether 
or not a design with nonover lapping RDCSSs exists is challenging [see, e.g., 
Bingham et al. (2008)] and actually finding such a design when it exists can 
be difficult. In the next section, we develop conditions for the existence of 
such designs. 
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4. Existence of RDCSS. In Section 4.1, results are given that focus on 
the maximum number of disjoint subspaces of equal size that are contained 
in the effect space V = PG{p — 1, 2). New results are developed in Section 4.2 
for the more general setting of unequal sized RDCSSs. This latter case is im- 
portant in multistage experiments where the number of units (nt) in a batch 
at stage i is not the same for all i = 1, . . . , to. The results are constructive 
and thus allow experimenters to find designs in practice. 

4.1. RDCSSs, spreads and disjoint subspaces. In most applications, the 
number of stages, to, of randomization is pre-specified by the experimenter. 
Thus, if a set, S, of disjoint subspaces can be obtained, with |<S| > to, an 
appropriate subset of S can be chosen to construct the RDCSSs. 

Definition 1. For 1 < t < p, a (t — l)-spread of the effect space V is a 
set, S, of (t — l)-dimensional subspaces of V, which partitions V. 

In a (i — l)-spread, S, every element of V is contained in exactly one 
of the (t — l)-dimensional subspaces of S. A (t — l)-spread is said to be 
nontrivial if t > 1. Given that a (t — l)-spread, S, of V exists, the size of S is 
|5| = (2 P — l)/(2* — 1). A necessary and sufficient condition for the existence 
of a (t — l)-spread is that t divides p [Andre (1954)]. As a result, if p is a prime 
number, there does not exist any nontrivial (t — l)-spread. Nevertheless, 
since we are interested in only m disjoint subspaces, we need only for the 
maximum number of disjoint {t — 1) -dimensional subspaces contained in V 
to be at least m. This is called a partial (t — l)-spread in finite projective 
geometry settings. 

Definition 2. A partial (t — l)-spread 5 of the effect space V is a set 
of (t — l)-dimensional projective subspaces of V that are pair-wise disjoint. 

Theorem 3 gives necessary and sufficient conditions for the existence of a 
set S of disjoint subspaces, each of size 2* — 1, for t < p. It also provides the 
size of the minimum overlap in the case where there are no two subspaces 
that are disjoint. The proof of the theorem can be deduced from that of the 
more general setup of Theorem 6. 

Theorem 3. Let V be a projective space PG(p— 1,2) and let S± and S2 
be two distinct (t — 1)- dimensional projective subspaces of V 0/ size \S~\\ = 
\S 2 \ = 2* - 1 for < t<p. Then: 

(a) If t < p/2, then subspaces S\ and S2 exist such that Si n S2 = 4>i 

(b) // t > p/2, then for any choice of S 1 ,S 2 in V, \S X D S 2 \ > 2 2t ~P - 1, 
and there exist subspaces Si, S 2 such that the equality holds. 
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Theorem 3(a) guarantees that, when t<p/2, one can obtain at least two 
disjoint (t — l)-dimensional subspaces of V . When p is not divisible by t, 
it can be expressed as p = kt + r for nonnegative integers k,t,r satisfying 
< r < t < p and k > 1. Since t divides kt, there must exist a (t — 1)- 
spread, S , of PG(kt - 1,2) contained in where |«S | = (2 fei - l)/(2* - 1). 
Then, since Sq is a (t — l)-spread of a subspace in "P, the set of disjoint 
(t — l)-subspaces in V may be expandable. The following result ensures the 
existence of a larger set of disjoint (t — l)-dimensional subspaces. 

Theorem 4 [Eisfeld and Storme (2000)]. Let V be a finite projective 
space PG(p — 1,2), with p = kt + r /or < r < t < p. Then, there exists a 
partial (t — l)-spread S of V with \S\ = 2 r \tZi — 2 r + 1. 

Next, Theorem 5 summarizes results available on the maximum number 
of pair-wise disjoint (i — l)-dimensional subspaces of V for different combi- 
nations of t and p. 

Theorem 5 [Govaerts (2005)]. Let V be a projective space PG(p-l, 2) , 
with p = kt + r for k > 1, < r < t < p, and let S be a partial (t — l)-spread 
of V with \S\ = 2 r \tZi ~ s > where s is known as the deficiency. Then: 

(a) I/r = l 5 then s > 2 r - 1 = 1; 

(b) // r > 1 andt> 2r, then s > 2 r ~ 1 - 1; 

(c) Ifr> \ andt< 2r, then s > 2 r ~ l - 2 2r ~ t ~ 1 + 1. 

Although Theorem 5 provides an upper bound for the number of disjoint 
subspaces in V, it does not guarantee their existence. As illustrated in the 
following two examples, the bounds may not be tight. 

Example 3. Consider a single-replicate 2 5 factorial experiment with 
randomization restrictions defined by 51,52 and S3 such that S\ D {A,B}, 
S2 D {C} and ^3 D {D,E}. As discussed in Section 3.2, for analyzing this 
design, we need at least three half-normal plots depending on the overlapping 
pattern among the S^s. The half-normal plot approach requires at least 
seven effects for each plot [Schoen (1999)]. Therefore, since the Si's are 
projective subspaces, the desired features of the RDCSSs are that Si CiSj = 4> 
for all i 7^ j and > 7 = 2 3 — 1 for all i. Since p = 5 is prime, there does not 
exist a nontrivial (t — l)-spread of V = PG(4, 2) with t > 3. However, p can 
be expressed as p = kt + r with t = 3, /c = l,r = 2 and, from Theorem 5(c), 
|<S| < 2. There is no certainty from Theorem 5 regarding the existence of 
even one pair of disjoint two-dimensional subspaces. In fact, from Theorem 
3(b), any pair of two-dimensional subspaces in PG(4,2) has an overlap of 
at least 2 2t ~ p — 1 = 1 effect; hence, the bound is not tight. 
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Example 4. Consider a 2 full factorial design with m stages of ran- 
domization characterized by RDCSSs defined by Si, ... , S m . From Theorem 
4, with p = 8, t = 3, r = 2 and k = 2,we know that there exists a partial two- 
spread where the number of disjoint two-dimensional subspaces of the effect 
space V is 33. Theorem 5(c) gives the size of maximal partial two-spread of 
V to be bounded above by 34. Thus, either the bound in Theorem 5(c) is 
again not tight or there exists a larger number of disjoint two-spaces of V 
than guaranteed by Theorem 4. 

Although the results presented in Theorems 3, 4 and 5 focus on PG(p — 
1,2), they can easily be generalized to PG(p — l,q), where q is a prime or 
prime power. See Eisfeld and Storme (2000) and Govaerts (2005) for the 
generalizations of Theorems 4 and 5, respectively. Theorem 3(a) also holds 
for two (t — l)-dimensional subspaces of PG(p— l,q). For the generalization 
of Theorem 3(b), if t > p/2 and Si, S2 are two (t — 1) -dimensional subspaces 
of PG(p - l,q), then |Si n S 2 | > (q 2t ~ p - l)/(q - 1) and there exist a pair 
such that the equality holds. 

We note that, for the special case of 2-level factors, t = 2 and p odd, we 
have p = 2k + 1 for positive integer k; so, r = 1 and Theorem 4 guarantees 
that the bound of Theorem 5(a) is achieved and |<S| = (2 P — 5)/3. For this 
case, a construction was proposed by Wu (1989) based on the existence of 
two permutations of the effect space satisfying certain properties. The result 
of Theorem 4 is more general since it holds for any integer t (0 < t < p) and 
is easily extendable for arbitrary prime or prime power q in PG(p — l,q). 

4.2. RDCSSs with subspaces of different size. So far, the results have 
focused on the existence of disjoint subspaces of the same size. However, it 
is not unusual for disjoint subspaces of different sizes to be required [see, 
e.g., Example 4 in Bingham et al. (2008) and the battery cells experiment 
in Vivacqua, Bisgaard and Steudel (2003)], and we next prove a new result 
that gives conditions for the existence of such a set of disjoint subspaces. 

Theorem 6. Let V be a projective space PG(p — l, 2) and Si be a (tj — 1)- 
dimensional subspace of V , where <ti <p for i = 1,2. Then: 

(a) If t\ + ti < p, then there exist S± and S2 such that Si D S2 = 4>; 

(b) Ifh+h>p, then for any choice of Si, S 2 inV, |SinS 2 | > 2 tl+t2 ~P - 
1, and there exist Si and S2 such that the equality holds. 

Proof, (a) Define the effect space V = {F±, F p ), where Fj's denote 
the main effects of the independent factors of a 2 P factorial experiment. Since 
h+t 2 <p, define Si = (Ft, . . .,F tl ) and S 2 = {F tl+ i, . . .,F tl+t2 ). Clearly, Si 
and S 2 are disjoint. 
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(b) For the case t\ + t 2 > p, define the subspaces Si and S 2 to be £1 = (F\, 
...,F tl ) and S 2 = (F p _t 2 +i, ■ ■ ■ , Ft i: F tl+ i, . . . , F p ). Thus, Si nS 2 = (F p _ t2+ i, 
...,F tl ) with |5in5 2 | = \PG(t 1 +t 2 -p-l,2)\ = 2 tl+ * 2 -P-l. Now, if there 
exist subspaces Sf and S 2 for which ti + t 2 >p and \S* n S 2 1 < 2 tl+t2 ~ p - 1, 
then it can be shown that | (5*1 , ^2) I > 2 P — 1 [see Ranjan (2007), Theorem 
4.2 for details]. This contradicts the fact that if S^ C V and S 2 C V, then 
(Sl,S 2 ) is also a subspace in V. Thus, S± and S 2 provide the minimum 
possible overlap, as required. □ 



For t\ = t 2 = t, Theorem 6 simplifies to Theorem 3. When t\ + t 2 < p 
[as in Theorem 6(a)], in addition to Si and S 2 , one can expect mutually 
disjoint subspaces of size 2* — 1 that do not overlap with Si and S 2 , where 
t <p — max(ti , t 2 ) . The next theorem, which is the main result of this section, 
establishes the existence of a set of unequal sized subspaces of V, where 
tj + tk <p for any pair of subspaces Sj and Sk ■ 

Theorem 7. Let V be a projective space PG{p — 1, 2) and Si be a (t% — 
1)- dimensional subspace of V . If p/2 <t\ < p, then there exist S 2 ,---,S m 
such that Si fl Sj = <f> for all i,j 6 {1, . . . ,m}, where \Si\ = 2** — 1 for ti < 
p — ti, 2 <i <m and m = 2 tl +1. 

Proof. Let p/2 <ti<p and define t* =p — t\. Then, the effect space 
V is a FG{ti + 1* - 1, 2). Let V' = PG(2ti - 1, 2) be such that V D V. Let 
Si be a (ti — l)-dimensional subspace of V, and let <S' be a (t% — l)-spread 
of V' that contains S\. Then, 

|5| - m - Ipg^-1,2)! -^3T" 2 +L 

The set of disjoint (t* — l)-spaces of V whose elements are disjoint from Si 
is given by <S* = {S fl V : S € S' \ {Si}}. The elements of <S* can be denoted 
by S 2 ,...,Sm- For every i = 2,...,m, since the desired subspace Si is of 
size 2 li — 1 with ti < p — ti = t* , one can construct Si by constructing a 
(ti — l)-dimensional subspace of S* . □ 

This theorem is important as it guarantees the existence of 2* 1 + 1 disjoint 
subspaces of different sizes, with one (t\ — 1) -dimensional subspace (p/2 < 
ti < p) and 2* 1 subspaces of dimension up to (t — 1) each, where t<p — t\. 
Note that the boundary conditions ti=p and ti =p/2 are uninteresting, 
as t\ =p implies that S\ is the entire effect space, and ti =p/2 guarantees 
the existence of a (ti — 1) -spread of V. The proof of Theorem 7 leads to 
a construction strategy for m = 2* 1 + 1 disjoint subspaces of unequal sizes, 
as illustrated in Section 5.3. Although Theorem 4 is not a special case of 
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Theorem 7, a similar construction is applicable for equal sized RDCSSs and 
is discussed in Section 5.2. 

Both Theorems 6 and 7 can be generalized to (/-level designs for prime 
or prime power q. Theorem 6(a) also holds for (ij — l)-dimensional sub- 
spaces Si's in PG(p — 1, q). For Theorem 6(b), if t\ + 1 2 > p, then \S\ fl S2I > 
(gti+*2-p — l)/(qr — l)s an d there exists a pair Si, $2 such that the equality 
holds. The generalization of Theorem 7 guarantees the existence of a set of 
m = q fl + 1 disjoint subspaces of different sizes, where |Sj| = (q li — l)/(q — 1) 
for p/2 < t\ <p and t{ <p — t\. The proofs of results for arbitrary prime or 
prime power q are similar to the q = 2 case shown above. 

Thus far, we have established some necessary and some sufficient condi- 
tions for the existence of a set of disjoint subspaces of the same size and also 
of different sizes. If the desired number of stages of randomization is less 
than or equal to the number of subspaces guaranteed to exist, we can obtain 
an appropriate subset satisfying the randomization restrictions required by 
the experimenter. Finding an actual design with the required properties is 
the next issue for the experimenter, and this is discussed in the next section. 

5. Construction of disjoint subspaces. A construction approach is now 
proposed for factorial designs with m levels of randomization. First, the 
construction for equal sized disjoint subspaces is presented, followed by the 
construction of disjoint subspaces of different sizes. The subspaces them- 
selves have no statistical meaning until the factors have been assigned to 
columns of the model matrix or, equivalently, to points in PG(p — 1,2). 
The set of disjoint subspaces obtained from an arbitrary assignment may 
not directly satisfy the experimenter's requirements for the factor levels. 
Consequently, we propose an algorithm that transforms a set of subspaces 
obtained from the RDCSS construction to a set of subspaces satisfying the 
desired restrictions on the factor levels. 

5.1. RDCSSs and (t — l)-spreads. When t divides p, we know that there 
exists a (i — l)-spread ofV = PG(p — 1,2) [Andre (1954)]. The construction 
of the spread starts with writing the 2 P — 1 nonzero elements of GF(2 P ) in 





Table 1 






Cycles of length N 




Si 


s 2 


Sn 





1 


TV — 1 


N 


N + l 


2/V-l 


(0-l)N 


(0-l)N + l 


ON -I 
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cycles of length iV [Hirschfeld (1998)]. An element w is called primitive if 
{w l : i = 0, 1, . . . ,u — 2} = GF(u) \ {0}. Let w be a root of a primitive polyno- 
mial of degree p for GF{2 P ). Then, the 2 P — 1 elements of the effect space V 
or, equivalently, the nonzero elements of GF(2 P ) are w l , i = 0, . . . , 2 P — 2. The 
element w l can be written as a linear combination of the basis monomials 
w°, ... , u; p . The element 

(7) w % = aow p ~ 1 + Qii(; p_2 + • — I- a p _2W + a p _i 

represents an r-f actor interaction 5 = (cto, a.\, . . . , a p _i) for Oj G GF{2) if ex- 
actly r elements in <5 are nonzero. Following this representation for the fac- 
torial effects in V, and using shorthand notation (iN + j) to denote w lN+ i 
for < i < — 1, < j < N — 1 (where is the number of cycles), cycles of 
length iV can be written as in Table 1. The columns define the projective sub- 
spaces Sx, . . . ,Sn. A necessary and sufficient condition [Hirschfeld (1998), 
Chapter 4] that there exists a (t — l)-space of cycles of length iV which is 
smaller than \PG(p — 1,2)| is that the greatest common divisor, gcd(t,p), of 
t and p is greater than one. Then, N = \PG(p — l,2)\/\PG(l — 1,2)|, where 
/ = gcd(t,p)). Thus, when t divides p, there exist 2* — 1 cycles each of length 
iV which lead to the formation of the required subspaces Sx , ■ ■ ■ , Sn of V 
with Si fl Sj = 4> for i ^ j; that is, Sx, ■ ■ ■ ,Sn constitute a (t — l)-spread S 
of the effect space V = PG(jp - 1, 2). 

A (t — l)-spread of PG{p — 1,2), obtained as above, distributes all the 
factor main effects evenly among the |<S| disjoint subspaces. However, re- 
strictions on the m stages of randomization are usually pre-specified by the 
experimenter and, as illustrated in Example 5, an RDCSS for a block design 
will contain no main effects whereas, for a split-lot design, several factor 
main effects may be assigned to one or more RDCSS. 

Example 5. Consider a single-replicate 2 6 factorial experiment with 
the randomization structure determined by a blocked split-lot design, where 
the experiment has to be performed in blocks of size eight each. Suppose 



Table 2 

The two-spread, S , of PG(5, 2) using the primitive polynomial w 6 + w + 1 
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s 3 
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EF 


DE 
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AEF 


DF 


CE 


BD 
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CDEF 
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ABCD 


ABCEF 


ABDF 
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BF 


AE 


DEF 


CDE 


BCD 
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ABEF 


ADF 


CF 


BE 


AD 
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BDE 


ACD 


BCEF 


ABDE 


ACDEF 
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ADEF 
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that the experimenter wishes to specify the factorial effects, ABC, BDE and 
CEF, to be confounded with the blocks (i.e., S? = {ABC, BDE, CEF)). In 
addition, suppose the trials proceed in a two-step process, where the restric- 
tions imposed by the experimenter on the two steps of randomization are 
such that S| ^> {A, B} and S| D {D}. As a result, there are three restrictions 
on the randomization of the experiment: one due to blocking the experimen- 
tal units (SI) and the other two due to splitting the experimental units into 
sub-lots (S2, S|). To use half-normal plots, it is desirable to have three dis- 
joint subspaces each of at least size seven (i.e., t = 3), where the subspaces 
should satisfy the restrictions defined by S%, S% and Sf. Here gcd(t,p) = 3, 
so N= \PG(5, 2)\/\PG(2, 2)| = (2 6 - l)/(2 3 - 1) = 9, and there exist 7 cycles 
of length 9 or, equivalently, 9 disjoint subspaces of size 7 (i.e., a two-spread 
of V). The two-spread S = {Si, . . . , S9} obtained using the primitive poly- 
nomial w e + w + 1, root w and representation (7) is shown in Table 2. Notice 
that each two-dimensional subspace of V in S contains at most one main 
effect, and none matches the design requirements. It is shown below how 
to transform the spread S to S* , such that S* contains the three disjoint 
subspaces S\, S| and S3* satisfying the experimenter's requirement. 

For the transformation of spreads, we use an appropriate collineation [see, 
e.g., Batten (1997)] of the projective space V. A collineation of PG(p — l,q) 
is a permutation / of its points such that (t — l)-dimensional subspaces 
are mapped to (t — l)-dimensional subspaces for 1 < t < p. The existence 
of a collineation / from S to S* is equivalent to the existence of a p x p 
matrix A4 such that, for every given S £ S, there is a unique S* € S* and, 
for every z € S, there exists a unique z* S S* that satisfies z*' = z'Ai. Note 
that the transformation of a spread amounts to relabelling the columns of 
the model matrix. As a result, one cannot find a collineation matrix A4 if 
the experimenter's requirement is not feasible. Moreover, if the desired set 
of subspaces is nonisomorphic to the spread we started with, then also there 
does not exist any relabelling to obtain the desired design. However, finding 
an appropriate collineation matrix, whenever it exists, is also nontrivial. 
Next, we propose an algorithm that finds a collineation matrix A4, if it 
exists, and concludes the nonexistence if one does not exist. 

The proposed algorithm is illustrated through the setup of the 2 6 exper- 
iment of Example 5. To obtain a set of disjoint subspaces satisfying the 
restrictions imposed on the three stages of randomization, we have to find 
an appropriate 6x6 collineation matrix A4. The proposed algorithm for 
finding the matrix M is outlined as follows: 

1. Select one of the (3) possible choices from the spread S in Table 2 to be 
a set of three disjoint subspaces. For example, suppose that Si, S3 and 
S7 are chosen for S|, S3 and Sj\ respectively. 
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2. Choose two effects from S\, one effect from S3 and three effects from £7 
to relabel these to the desired effects (A, B),D and (ABC, BDE, CEF) in 
52,5*3 and S\, respectively. For example, one choice among the ( 2 )(i)( 3 ) 
different options is {CDE, BCF , D, EF , AC , BF}. The collineation ma- 
trix is defined by the mapping induced by CDE — > A, BCF — > B, D — > 
D,...,BF-> CEF. 

3. Denote the (i, j)th entry of the pxp matrix A4 as Xk for k = j + (i — l)p 
(i.e., list the elements row by row). Construct a p 2 x p 2 matrix Q and a 
p 2 x 1 vector 5 as follows. Define the rows of matrix Q and vector 5 in 
the order of restrictions on the transformation. For the example under 
consideration, the first transformation, (CDE) 1 AA = (A)', can be written 
as 

(8) [00111 0]M = [1 0]. 

In total, there are p (1 < s < p] independent transformations, and for 
each transformation there are p (p(s — 1) < i <ps) rows of 5 and Q. The 
first set (s = 1) of p = 6 elements of 5 are given by the right-hand side of 
(8). The corresponding rows of Q are defined as 

Qu = l, if / = (r-l)p + mod(i-l,p) + l 

and the rth entry of [0 1 1 1 0] is nonzero, 

= 0, otherwise, 

for 1 < i < p. Similarly, all the rows of the matrix Q and the vector 6 can 
be expressed using p restrictions on the transformation. For example, the 
second set (s = 2) of p = 6 rows (5i, Qn for p+ 1 < i < 2p) are given by the 
right-hand side of (BCF)'M = (B)'; that is, [0 1 1 1]M = 
[0 1 0], and so on for the remaining transformations. 

4. If there exists a solution for Qx = 5, reconstruct the pxp matrix M. (row 
by row) from the p 2 x 1 solution vector x = Q~ L 5, where Q~ L is a left 
inverse of Q. Exit the algorithm. 

5. If there does not exist a solution then go to step 2 and, if possible, choose 
a different set of effects from the subspaces selected in step 1. 

6. If all possible choices for the set of effects from these three subspaces 
have been exhausted, go to step 1 and choose a different set of three 
subspaces. 

7. If all the (3) choices for a set of effects have been used and still a solution 
does not exist, then either the two spreads S and S* are nonisomorphic 
or the experimenter's requirement is not achievable. Thus, the desired 
spread cannot be obtained from S. 

In the illustration used above, the effects chosen for relabelling the columns 
to achieve the desired design do provide a feasible solution to Qx = 5. The 
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collineation matrix A4, obtained from the solution x = Q L S, is given by 

1\ 


1 

o • 
1 

1/ 

In this example, an exhaustive search found that 45.7% of all possible 
choices give a feasible solution to the equation Qx = 5. That is, an arbitrary 
choice of p independent effects from S (according to steps 1 and 2) results 
in a feasible design 45.7% of the time. The rest of the time, the solution is 
infeasible since the full factorial design becomes a replicated fraction. 

The spread acts as a template which enables a faster search than the ex- 
haustive relabelling of all the factorial effects to find the design satisfying 
the experimenter's requirement. For this example, our algorithm requires 
at most (g) (g) ([) (3) ~ 5 x 10 5 different relabellings, whereas an exhaustive 
relabelling approach can require up to (2 6 — 1)! ~ 2 x 10 87 different rela- 
bellings. While searching through each of the relabellings in the template 
can be time consuming, our Matlab 7.0.4 implementation of the algorithm 
found the first feasible collineation matrix in 5.34 seconds on a Pentium(R) 
4 processor machine running Windows XP. 

5.2. Partial (t — l)-spreads. When t does not divide p, Theorem 4 guar- 
antees the existence of \S\ = 2 r \tZ\ — 2 r + 1 disjoint subspaces, where 
p = kt + r. For constructing these subspaces, one can use the existence proof 
of Eisfeld and Storme (2000) for the most part. However, the proof as- 
sumes the existence of an (sj — l)-spread S[ of V% that contains an (s.; — 1)- 
dimensional subspace, C/j, of V i+ i, where Sj = it + r, Vi = PG(2si — 1,2) 
and V i+ i = PG(si + t — 1,2) for i = 1, . . . ,k — 1. This is nontrivial, and we 
propose a two-step construction method to get this as follows: (a) construct 
an (si — l)-spread S" of Vi as described in Section 5.1 and then (b) trans- 
form the spread S" to S[ by finding an appropriate collineation such that 
Ui £ 5,-. Thus, we can construct a set of |5| disjoint (t — l)-dimensional 
subspaces. Finally, we find an appropriate collineation to transform the par- 
tial (t — l)-spread S to obtain the m RDCSSs satisfying the experimenter's 
requirement. 

5.3. Disjoint subspaces of different sizes. A more general setting is when 
the RDCSSs are allowed to have different sizes. For a 2 P full factorial design, 
Theorem 7 guarantees the existence of one subspace S\ of size 2* 1 — 1 
with t\ >p/2 and 2* 1 subspaces of size at most 2 p ~ tl — 1. For constructing 
these 2* 1 + 1 mutually disjoint subspaces of V, the proof of Theorem 7 
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Table 3 

Three-spread of PG{7,2) containing Si = (A,B,C,D), S 2 D {E,F} and S 3 D {G} 
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requires constructing a (t\ — l)-spread S" of PG(2t\ — 1,2) that contains 

51. The spread S" can be obtained by first constructing a (ti — l)-spread of 
PG(2t\ — 1, 2) and then by applying the appropriate collineation {Mq) found 
by the algorithm described in Section 5.1 to obtain S 1 that contains Si. After 
S = {S nV : S € S' ' \ {Si }} is obtained, one has to find a suitable collineation 
(Mi) so that the final set of subspaces satisfy the experimenter's restrictions 
on RDCSSs. The algorithm can be made more efficient by combining the 
problem of finding the two collineation matrices into one. When transforming 
the spread S" to S' containing Si, we can impose other restrictions on 

52, --.,S m in this step itself. The steps of the construction are illustrated 
through Example 6. 

Example 6. Consider a single-replicate 2 7 factorial design with three 
stages of randomization. Let the restrictions imposed by the experimenter 
on the three RDCSSs be Si D {A,B,C,D}, S 2 D {E, F} and 5 3 D {G}. 
Following the notation of Theorem 7, since p = 7 and t± = 4, there exist m = 
2 4 + 1 = 17 pair-wise disjoint subspaces with \Si \ = 2 tl — 1 for i = 1, . . . , 17, 
where ti = 4 and ti < 3 for i = 2, . . . , 17. First, we construct a three-spread 
S" of PG(7,2) using the method described in Section 5.1, and then we find 
an appropriate collineation matrix Mq that transforms S" to S such that 
S contains Si = (A,B,C,D), S2 D {E,F} and S3 D {G}. Table 3 contains 
some of the elements of the transformed spread. 

For this example, the three disjoint subspaces Si,S2 and S3 that satisfy 
the experimenter's requirements are obtained by deleting the elements of 
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Si's in the transformed spread that contains H. The collineation matrix Ai 
used for the transformation is as follows: 
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Thus, {S nV:S € S' \ {Si}} U Si contains the required set of subspaces 
S\,S2 and S3 for the three stages of randomization. 

6. Fractional factorial designs. Though the designs discussed in this ar- 
ticle are full factorial designs, the results also apply to FF designs with 
randomization restrictions. The existence of regular FF designs with ran- 
domization restrictions is equivalent to that of the full factorial design gen- 
erated from the basic factors of the FF design (base factorial design). One 
can find a 2 r ~ s fractional factorial with randomization restrictions by first 
investigating the existence a 2 U full factorial design (u = r — s) with the ap- 
propriate randomization restrictions. If the design exists, one can construct 
the 2 U full factorial design with randomization restrictions as outlined in 
Section 5 and then assign the added factors to the effects at the appropriate 
level of randomization to get the desired design. 

Consider, for example, a 2 8-2 fractional factorial experiment with ran- 
domization structure characterized by a split-lot design. Suppose that the 
experiment has to be run in four stages with randomization restrictions 
given by Si D {A,B}, S 2 D {C,D}, S 3 D {E, F} and S 4 D {G,H}. Then, 
the six independent basic factors (A,B,C,D,E,F) with effect space V = 
(A,B,...,F) form a 2 6 full factorial split-lot design. The results discussed 
in Section 4 guarantee the existence of a two-spread of V, and the algo- 
rithm described in Section 5.1 shows the construction for obtaining three 
disjoint subspaces of size seven each satisfying the requirements of Si, S2 
and S3. Moreover, we know that there exist nine disjoint subspaces of size 
seven each. Therefore, S4 can be constructed by choosing a subspace from 
the remaining six disjoint subspaces and then by aliasing two generators G 
and H with effects in this subspace [Box, Hunter and Hunter (1978)]. 

Fractionation of the base factorial design can occur in a different scenario. 
For example, consider a 2 8-2 experiment with the requirement of three stages 
of randomization and the restrictions on the RDCSSs be defined by Si D 
{A,B}, S2 D {C,D,E} and S3 D {F, G, H}. In this case, one can use the 
algorithm in Section 5.1 to construct three disjoint subspaces satisfying Si D 
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{A,B}, S2 D {C,D,E} and S3 D {i 7 }- Next, one can choose two generators 
(or points) from S3 that are assigned to G and H . 

To rank the designs, one can use existing criteria such as minimum aber- 
ration [Fries and Hunter (1980)], maximum number of clear effects [Chen, 
Sun and Wu (1993) and Wu and Chen (1992)] and ^-criterion [Bingham et 
al. (2008)]. Thus, one can select an appropriate set of generators based on 
the experimenter's interest. The structure of a spread acts as a template to 
shorten the computer search for good fractional factorial designs. 

7. Concluding remarks. In this paper, we have demonstrated that the 
projective subspaces of the effect space V can be used to characterize the 
randomization restrictions of factorial designs in block, split-plot, split-lot 
and other structures. Under the assumptions of model (1), Theorem 2 sum- 
marizes the impact of randomization restrictions on the distribution of fac- 
torial effect estimators and motivates the search for disjoint randomization 
defining contrast subspaces. Obtaining a set of disjoint subspaces of the ef- 
fect space V is nontrivial. In the most general case, Theorem 7 presents a 
necessary and sufficient condition for the existence of disjoint subspaces of 
unequal sizes. Furthermore, the proof motivates strategies for constructing 
these, and this makes a wide variety of designs accessible in this unified 
framework. 
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